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ABSTRACT 


A Nusselt type analysis was performed during laminar film condensa- 
tion on the inside of a rotating truncated cone. This analysis was 
employed to determine the condensing limit of a wickless heat pipe, 
rotating about its longitudinal axis. 

Performance characteristics including the effects of geometry, 
rotational speed, and the characteristics of fluid are given. A 
comparison is made between the condensing, boiling, sonic, and entrain- 


ment limits for a given heat pipe geometry. 
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I. INTRODUCTION 


A. BACKGROUND 

In recent years considerable interest has been generated in the 
study of capillary-pumped heat pipes as a means of transferring large 
quantities of heat. At present, for most pipes using ordinary fluids, 
the limit on the heat transfer appears to be the rate that the condensate 
can be returned to the evaporator [1]. The capillary-pumped heat pipes 
are primarily limited by the small pressure differentials obtainable 
across the meniscus of the liquid in the wicking, which gives rise to 
the pumping force. This differential pressure must overcome the losses 
of the liquid flowing in the wick, the losses in the vapor flow, and, 
depending on the orientation, the gravitational forces. A second limit, 
the onset Of nucleate boiling in the evaporator portion of the wicking 
structure, may disrupt the capillary flow [2], [3]. Nucleate boiling will 
cause the formation of bubbles in the liquid condensate. These bubbles 
may then become trapped in the wicking structure, thus impeding the re- 
turn flow and limiting the operation of the pipe [2], [3]. 

The limitations of the capillary~pumped heat pipes can be overcome 
through the use of centrifugal acceleration [4]. A pipe, rotating about 
its longitudinal axis, will generate its own gravitational field, indepen- 
dent of orientation, and thereby force the condensate to return to the 
boiler. 

A rotating heat pipe is a closed hollow shaft which contains a work- 
ing fluid. Rotation about its longitudinal axis generates a rotational 
force field. The rotational acceleration, together with a slight internal 


wall taper, pumps the working fluid from the condenser to the boiler 


section. There is no porous wick structure as in a capillary heat 
pipe. The operation of a rotating heat pipe is shown schematically 


im Fig. le 
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Figure |* Cutaway View of Typical Rotating 
Heat Pipe 


As the pipe is rotated, an annulus of working fluid is formed at 
the largest internal radius of the shaft. Heat is applied to this 
portion of the pipe. As heat is applied, the fluid evaporates or 
boils, forming a vapor. This vapor is then at a higher pressure than 
that at the narrower end of the pipe (condenser). This pressure differ- 
ential forces the vapor from the boiler to the condenser section. Heat 
is removed at the condenser section, causing the vapor present at that 
section to condense on the walls of the shaft. The centrifugal force 
generated by the rotational motion of the pipe insures that the con- 


densate will return to the boiler section. Thus the rotating heat 
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pipe has a closed, two-phase heat transfer cycle, and the flow of the 
condensate is not retarded by the porous wick structure. 

The rotating heat pipe may be applied wherever there is a rotatirg 
part and a need to transfer heat. A few applications of the rotating 
heat pipe are the cooling of: a high speed drill, the rollers for a 


conveyer table of hot metals, the rotors of electric motors, and turbine 


blading. 


B. THESIS OBJECTIVES 
The objectives of this study were: a) to understand the operational 
principles, b) to develop an analytical model and, c) to predict the 


performance characteristics of a rotating heat pipe. 


Ll 


II. LIMITS OF OPERATION 


As in the capillary heat pipes, the heat transfer limits of the 
rotating heat pipe are controlled primarily by fluid dynamic considera- 
tions.[{[5]. The effect of the following limits will be discussed: a) 
the boiling limit, b) the entrainment limit, c) the sonic limit, and 


d) the condensing limit. 


A. BOILING LIMIT 

When a surface is heated in a stagnant pool of saturated liquid, 
three basic regimes of heat transfer are observed to occur. In the low 
range of temperature differentials ,(T,-T,) ,the heat flux is low and the 
principle mechanism of heat transfer is free convection. As the tempera- 
ture differential is increased, nucleate boiling becomes dominant. In 
this region, bubbles of the vapor are formed on the heating surface, 
grow in size, and eventually break off. After breaking away they rise 
to the surface and burst, releasing the vapor. The formation of the 
bubbles and their consequent rising agitates the liquid, causing convec- 
tive currents to be set up. This agitation intensifies the heat transfer 
mechanism causing an upturn in the heat flux curve and increasing the 
heat transfer coefficient significantly. As the temperature differential 
is increased further, the bubbles become so numerous, that they begin to 
coalesce, and a film of vapor can then blanket the surface. The insu- 
lating effect of this vapor blanket makes the film boiling regime un- 
desirable. A maximum heat flux exists for every boiler surface at which 


the transition occurs between nucleate boiling, and film boiling. 
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During the operation of the rotating heat pipe, the heat transfer 
coefficient, as mentioned above, will be increased due to nucleate boiling. 
A further increase in the coefficient will be experienced due to the ro- 
tational acceleration of the boiler [6], [7], [8]. The boiling limit in 
the rotating heat pipe will be reached upon the departure from nucleate 
boiling. This will occur as a sudden drop in the heat transfer coef- 
ficient which occurs when the maximum heat flux mentioned above has been 
reached. If the power input is held constant during this time, the sur- 
face temperature of the container will rapidly rise, usually causing the 
surface to melt. In the literature this phenomenon has been termed 
departure from nucleate boiling, boiling crisis, or burnout [9]. The 
boiling limit of the rotating heat pipe may be derived directly from 
the work done on conventional non-rotating boilers. 

Kutateladze has postulated that the transition from nucleate boiling 
to film boiling is purely a hydrodynamic phenomenon [10]. Based on this 
assumption he has obtained the correlation formula for predicting the 


burnout heat flux. 


V4 
=K Jp, Mal 2 al P, -p |} aa 


V 


pe 


| 


> 


b 


where 

q = total heat transfer, BTU/hr 

A. = heat transfer area in the boiler, a 

fa = density of the vapor, tbat 

h = latent heat of vaporization, BTU/1bm 
= surface tension, lbf/ft 

‘ = acceleration of gravity, cee 


P , - density of the fluid, ene 


Kutateladze's experimental data indicates the average value of K is .14, 
obtaining from a range of .13 to .19 for various surface conditions of 
horizontal wires and discs [11]. A counterflow configurattion of in- 
viscid vapor and liquid jets in the high nucleate boiling regime was 
proposed by Zuber, Tribus, and Westwater to predict the burnout point [12]. 
Their results for the burnout heat flux on a flat plate heater shows that 
K is equal to 77/24 (or.i3). Experimenters have substantiated Kutateladze's 
formula and have shown that the burnout heat flux increases for rotation 
approximately as the rotational acceleration to the one fourth power [13], 
[14]. Thus the boiling limit of the rotating heat pipe may be adequately 
described by the correlation formulas of Kutateladze and Zuber [9], [15]. 
The heat flux equation (1) may be recast as a heat rate by multiplying 
through by the surface area of the rotating heat pipes boiler. This gives 


equation (2) 
1/4 
9, = K vp, be aoa | A -py] } (2) 


where 


Bi ENTRAINMENT LIMIT 

During the operation of the rotating heat pipe, a counter current 
flow of the liquid condensate along the walls and vapor in the central 
core of the pipe is set -. A shearing off of the condensate by the 
vapor,holding up the condensate and thereby causing the vapor formation 
in the boiler to exceed the condensate return, will be defined as the 
entrainment limit. High heat fluxes during steady state operation and 
certain critical heat fluxes during start up conditions may possibly 


generate counter current (vapor-condensate) flow rates capable of 


triggering entrainment of the condensate. In counter current flow, a 
condition closely related to the entrainment problem has been defined in 
the literature as flooding [16]. When this condition is reached, it is 
impossible to increase the flow rate of the liquid component without 
making a corresponding decrease in the flow rate of the vapor component 
and vice versa. Hewitt, Lacey, and Nicholls have shown that the flow 
contact length:‘is a strong parameter for the flooding transition at 
high liquid flow rates. They also have shown experimentally that as the 
absolute pressure is increased, the flow rates attainable before the 
flooding transition are higher [17]. In order to understand the state- 
ments and results of Hewitt, Lacey, and Nicholls, this discussion of 
droplet entrainment from a liquid film for co-current annular flow by 
C. B. Wallis is presented: 

"For any given liquid flow rate the gas-liquid interface is 
characterized by the following physical appearance. At zero gas 
velocity the liquid has a smooth, mirror-like appearance as it 
moves along the duct. At a low gas velocity the surface starts 
to break up into a ripple pattern. As the gas velocity is 
further increased, the ripples start to bunch up and a wave is 
produced which appears to move at a velocity several times that 
of the entrained film. This wave, referred to as a roll wave, 
moves down the duct on top of the entrained film. At still 
higher gas velocities these roll waves grow smaller (yet more 
frequent) and droplet entrainment becomes evident. These waves 
are a precondition for the onset of entrainment and the source 
of the droplets. The wave tops are sheared off by the gas 
flowing over them. The roll waves appear to ‘roll’ in the 
direction of the gas flow." [18] 

The rotational acceleration in the rotating heat pipe will be a 
significant parameter in suppressing the formation of a rolling wave as 
described by Wallis [18]. It should act in the same manner as the in- 
creased absolute pressure delayed the onset of flooding in the experiments 


of Hewitt, Lacey, and Nichol]s [17]. In addition, since the rotating heat 


pipe is a sealed, self contained unit, the internal pressure may be 


regulated to such a value as to ensure no liquid entrainment will 
occur. 

At the present time there is no satisfactory theory available to 
predict the onset of entrainment, but the experimental data of Hewitt, 


et al, can be used to predict the entrainment limit [17]. 


Cy SONIC LIMIT 

As the heat flux in the rotating heat pipe is increased, it is 
possible that the flow rate of the vapor down the container may be 
limited by the attainment of a choked flow condition for the vapor. 
This choking of the vapor flow will definitely limit the amount of 
energy the vapor can transport down the pipe, thus limiting the opera- 
tion of the pipe. The heat transfer mechanism employed in the rotating 
heat pipe depends on the latent energy of the working fluid. The heat 


transfer rate in the pipe may then be expressed as 
q, = A Ng (3) 


where 
my, = mass flow rate of the vapor (lbm/hr) 


From continuity, the mass flow rate of the vapor may be expressed by 
mM) = pyUA (4) 


where 


velocity of the vapor (ft/sec), and 


c 
li 


A = cross sectional area for flow (HE) 
Combining equations (3) and (4) yields 


= h. 
9 Pty ‘ 1Q (5) 
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If the cross sectional area of the pipe is known where the choked condi- 
tion ees. the sonic limitation on the heat rate may then be com- 
puted using equation (5) with the vapor velocity taken as the sonic 


velocity. 


D. CONDENSING LIMIT 

The condenser may limit the rate of return of the condensate to the 
boiler by two mechanism; a) the rate that energy may be extracted from 
the system, and b) the maximum flow rate of the condensate. At the 
present no solution exists in the literature for condensation on a 
truncated cone. The solution for condensation on a rotating cone has 
been carried out by Sparrow and Hartnett [19]. Their solution employs 
a boundary layer similarity approach. As is stated by Sparrow and 
Hartnett, the solution applies only to cones that are not too slender 
and in a region sufficiently removed from the apex of the cone. As can 
be seen in Fig. l, the rotating heat pipe will generally employ a 
slender, truncated cone for a condenser section. If Sparrow and Hartnett's 
solution for condensation on a cone were employed, then the boundary condi- 
tion of the rotating heat pipe would not be fulfilled, namely that the 
velocity of the condensate is zero at the narrowest section of the con- 
denser. Hence, a truncated cone solution is necessary. However, since 
a similarity solution for the truncated cone would not apply in the 
region near the narrow end of the condenser, due to the boundary layer 


assumptions, the following Nusselt-type analysis was made. 


1 

As slender pipes are anticipated (small cone angles) then as a 
good first approximation, it can be assumed that choked flow occurs at 
the exit area of the boiler. 


III. THEORETICAL ANALYSIS OF FILM CONDENSATION 
ON THE INSIDE OF A ROTATING-TRUNCATED CONE 
In attacking the problem of condensation on the inside of a trun- 


cated cone the following assumptions are made. 


A. ASSUMPTIONS 

a) Film condensation, not dropwise condensation, occurs in the 
condenser section. Thus the condensing liquid wets the container sur- 
face and condensation occurs at a sufficient rate to allow a continuous 
film over the surface. The assumption of filmwise condensation is felt 
to be the conservative choice as dropwise condensation has been shown 
to give higher heat transfer rates [20]. 

b) The film of condensate undergoes laminar flow. Cannon and 
Kays, using a pipe rotating about its longitudinal axis with annular 
flow, have reported, 

"Tube rotation apparently has the effect of stabilizing a 
laminar flow so that a higher through-flow Reynolds number is 
obtained before transition occurs." [21]. 

c) The fluid properties are constant. Heat pipes are normally 
operated as nearly isothermal devices. This condition should be easily 
achieved in the rotating system and thus the fluid properties will very 
nearly be a constant. 

d) The subcooling of the condensate may be neglected. A small 
correction factor is normally applied to the latent heat of vaporization 
term in order to take into account the effect of subcooling. The neces- 
sity for this correction is determined by the comparison of the term 
C aa ce is the specific heat at constant pressure,BTU/lbm deg F, and 


AT is the difference between the wall temperature and the saturation 
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temperature, deg F) with the latent heat of vaporization. Normally 


os is much less than h. and thus may be neglected. 


fg 

e) The momentum changes through the condensate are small and 
thus there is essentially a static balance of forces. Sparrow and Gregg 
have shown that the momentum effects are negligible for ordinary fluids 
in laminar film condensation on a flat plate [22]. Assumption i), a thin 
film, also indicates that momentum changes may be neglected. 

f) The vapor exerts no drag on the motion of the condensate; 
there is no interfacial shear. Experiments have shown that for Prandtl 
numbers of ten and higher the effect on heat transfer of interfacial 
shear is negligible. For Prandtl numbers near one the effect is only 
slight (within five eae (23). 

g) There is a linear distribution of temperature through the 
condensate film. Rohsenow has shown that for Prandtl numbers above one, 
the non-linear distribution of temperature through a thin film, due to 
energy convection, is hardly distinguishable [24]. Thus energy convec- 
tion is neglected and heat transfer occurs by conduction only. 

h) The centrifugal force is much greater than the force of 
gravity and thus gravity may be neglected. 

i) The thickness of the film is much less than the radius of 
curvature of the condenser wall. 

j) The vapor space is essentially at one pressure, Pie 
There is no radial distribution of pressure in the vapor space and the 
axial distribution of pressure is nearly a constant [25], [26]. 

k) The density of the fluid is much greater than the density of 
the vapor. 

1) The rotating heat pipe is operating at a steady state condi- 


tion. 


m) Velocity gradients in the circumferential direction rela- 


tive to the pipe wall are negligible. 


By using the above assumptions and the coordinate system shown in 


Fig. 2, an analysis similar to Nusselt's original film condensation 


theory may be used [27]. 





Condenser 


Figure 2° Coordinate System ond Geometry 


for Condensing Limit 


A force balance may be taken on an infinitesimal fluid element 


in the X-direction as shown in Fig. 3. This leads to an expression 


for the momentum equation within the condensate. 
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Figure 3 Force Balancein X-Direction on an 


Infinitesimal Fluid Element 


B. MOMENTUM EQUATION (X-direction) 


- oP 
dF, =2mr Pay —2nr {P+ ox} dy r2~r de + (6) 
OT is 2 2 
{r fi; i thee cds 5 2 Py wW dmwr ay ary. sin F 
where 

2 

P = pressure, (lbf/ft §) 

r = radius, (ft) 

X ™ co-ordinate measuring distance along the surface, (ft) 


Y = co-ordinate measuring distance normal to surface, (ft) 


tT = shear stress, (Ibé/£t7) 
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od = half cone angle (rad) 
Pe ™ density of fluid, (1bm/£t*), and 
WwW = angular velocity, (1/hr) 


This equation reduces to 


Gees Of 
O+-5% ay TPE! sing (7) 


C. MOMENTUM EQUATION (Y-direction) 
In like manner, using Fig. 4,a force balance in the Y-direction, 


neglecting the centrifugal head due to the vapor since py ee, » gives: 
p 

DF = 2arPdx - OTs P+ -~ dy} dx — pw fer cost dx dy (8) 

This reduces to 


= - = pure COS ‘f (9) 
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Figure 4° Ferce Balance in Y-Direction on an 


infinitesimal Fluid Element 


bad 


Making use of Fig. 2, the radius at any point in the fluid 


condensate may be expressed as 
r(xy)=Ro+X sing — Ycos f (10) 


where 
RS = minimum wall radius in the condenser section (ft) 


Substituting equation (10) into (9) for r(x,y) and integrating gives 


~N 


A é 
Vv 
dP 2 
fay 9-H cosh JF + Xsing ~ ¥ cose dy (11) 


where 


P = pressure of the vapor GbE Et) 


Ct 


I 


film thickness (ft) 


This yields 


P= P+ p,«icos ¢ {[Ro+ X sin | [SY] i . cosd 5 - ¥}} (12) 


Now differentiating equation (12) with respect to X gives 


6P dP. 7 WHS. d8 
a : = + pw cosd tRo- = 4 5sing ee ae (13) 


- Ysing -Scos pi} 


P 
If the result of equation c3), @e , is substituted into equation 


Ox 
(7) it becomes 
naoex wale ¢ | dd | ds : 3 $s) 
O = ae pics {Row Fasing +Xsinge Ysin «p Scosif ye (14) 
. , A 
t h bt j {R.+ f— LC lo = — 


Integrating equation (14) over the limits from Y to 0 and T to Oe 


then rearranging yields 





! di ~ ; c aes “" < 
er am bale “a oo | ee a (15) 
-Scos¢ +} {8-v} + pwising Ro dsing Yi 
where in equation (15) for laminar flow 
=a 
r= ie ie 
and 
u. = velocity of the fluid, (ft/hr) 
Be = viscosity of the fluid, (lbm/ft hr) 
Equation (15) may now be solved for the velocity profile ue 5 Sh 
integration over the limits 0 to Y and QO to Ue 
< GOSsp iS dd 
__ 1B pie 0) pcs foe 
Do ede ae, ji, bay Ssingds Xsin ip 5 (16) 


i Wwe ‘; w sing Y 
i Y>- 3-} a aaa {Ro+ xsing} {sy-} 


—O0%5 is 


Appiying assumption jj), ( vapor space at one pressure, PY ), allows 
equation (16) to be rewritten as 
P Ww = y) ( d ’ ‘ 18 
ee {AY ie a ml aos") 28) 4 xsin-flrn Pos) _ a?) 
j=, 4 ‘ Lag r 








Note that in equation OT aes ufx,y) since = Cx is 
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D. CONTINUITY EQUATION 


The continuity equation of fluid dynamicg requires that 


m= pp uA 
where 
m = mass flow rate of fluid, (lbm/hr) 
Ul = average velocity of fluid, (ft/sec) 


This can be written as 


(18) 
5 


Mg Jpcugoare ay 


a) 


Substituting equation (17) into (18) and carrying out the integration 


ee 
ys HEE rang Berea] S [Ryexsinf=Beossf} as 
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J 
© 4 


CTS 


lov» 


5 
FE rrcral- 5, 
Utilizing assumption i) (thin film) allows cquation (19) to become 


re a O ff. \ ; 4 i 
fo flg = {[Rot Xsin a oS) 


Differentiating this result with respect to X gives 


- 3 


asl 2 2 | 
J sin ho Ae te “3 SIN p SIR +X sind I} (21) 
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Ee ENERGY EQUATION 
If the film surface temperature is T. and the wall is at Tees then 


the heat transfer by conduction to an element of surface area dA O is 


kec OS p UE -T i dA. (22) 


dja= ——— 
44 {R rae ph { Rot X sing 
Sin Pp n : < 
0 
Rot Xsing — Scos¢4. 
where 
dA = C1) Rg t XSI fur tx 
A. = heat transfer area in the condenser, (ft) 
ke = thermal conductivity of the condensate film, 
(BTU/hr ft deg F) 
T. = saturation temperature, (deg F) 
T = inside wall temperature, (deg F) 
Therefore 
[> ae 
Oom«scos” | w_!A 
2 ie 3 j 23 
ng | { Rn tA Sal oe 
He a a \ : oa ‘ — 
Rot Asin arog fie 
which can be written as 
6 Cl I bl 1» 
| : | (24) 


J i ee aes 
Pn \| Rt xX nia s 


The assumption of a thin film, assumption i), gives 


Sere 
= — | ic i\| 
Roa a sin | 
Therefore 
: Sas] 5 coe 
Inf{y— aE SOR Oe 
Rit X sin] R id sin’ 
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Equation (24) then becomes 


daz gf 2m {Rr Xi j ti T- Tg a = 


The rate of condensation is 


dq,=h, dim, (26) 


fe 
Since the pipe is operating at a steady state condition, continuity 


for the entire system requires that 


and applying this result to equation (26) gives 
day = hy im, ey 


Equating equations (27) and (25), the mass rate of condensation per 


unit length is 


din KK? Cz t 
—f= = a r tute UG eh y (28) 
d x Shit, q L% Pj = fd) 
Equating equations (28) and (21), then rearranging gives 
a iy: * fBotxsind) i 


Bi iam Ps 


Be FILM THICKNESS 
New carrying out this integration (see Appendix 4) for 5 » the 


final results are 


zs Bix 4 
ie ) ce “ee? 





The local heat transfer coefficient, h, (BTU/hr ft 
by 


deg F), is given 





(31) 
er 


MEAN HEAT TRANSFER COEFFICIENT 


The mean heat transfer coefficient over the whole surface area 
of the truncated cone is given by 


Xe 

| ie ) i 
h= : wa eihia “a Ak 
m TLI2R +L sind ? ay 


(32) 


where 


Coa 


length along the wall of the condenser, (ft) 
Integrating equation (32) (see Appendix B) yields 


| 


: {< 
"Im {2k +L sin gp} 





(33) 


= V4 
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I¥. DETERMINATION OF THE CONDENSING LIMIT 


The analysis of film condensation on a truncated cone may be ap- 
plied to the condenser section of the rotating heat pipe. This then 
will give an equation which will predict the condensing limit for the 
rotating heat pipe. The total heat rate in the condenser, qe is given 
by 

ay = Brat Le G~ Tyh {2 Rot Laing} (34) 
The value of ho is now substituted into equation (34) and 


93 a3 
{{Ro+ L sin t) | ania § he (35) 





he Rare <a 
; 2 k¢ Pf w h 
aye {1 


This then is the operating equation for predicting the limiting 
heat rate for the condenser. In the equation, however, the quantity 
T, is not a readily controllable quantity. The inside wall temperature 
will actually adjust to a temperature corresponding to the heat rate 
(the saturation temperature) and the outside wall temperature. As the 
heat rate in the rotating heat pipe is increased, the internal pres- 
sure and saturation temperature will increase, accordingly requiring 
that the inside wall temperature adjust for this increase. Thus, the 
outside wall temperature of the condenser will be the controllable and 
measurable quantity. The inside wall temperature, Ts then may be found 
as a function of the outside wall temperature, Eos and the heat rate 


by use of Fourier's conduction law through the consenser wall, 


de, = ‘w2 miR+Xein ff Fete ee (36) 


where 


rk 


thermal conductivity of the wall, (BTU/hr ft deg F) 


*t = condenser wall thickness, ft 


Once again the cylindrical effects have been neglected due to the thin 
wall of thickness t. Integrating equation (36) from 0 to a and 0 to L. 


yields 


qi 2 LT, -T of oe sing} (37) 


A. OPERATING EQUATIONS 
The operating equation for the condensing limit is then an itera- 


tive solution between equation (35) and (37). 


3 it i i as a3, 
nee C aaa Pm ane 
aia um, 5 ay , LER, + Loin g| Ry! (35) 
k 
q, = mec nae -T HR, + 5 shee £ sing} (37) 


At this point it should be noted that although operating equa- 
tions for the condensing limit of the pipe itself have been derived, 
as an engineering consideration, the outside cooling mechanism for 
the pipe is extremely important. Indeed, it may become the limiting 


consideration. This then gives a third equation to be pat bel wal, 


= 
Sle ale LL, ala ot Dr sin: N (38) 


-this is assuming that the film and wall thicknesses are very thin, 


and thus the areas are approximately equal. 
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where 
ho = coolant heat transfer coefficient (BTU/hr na deg F) 


T, = coolant temperature, (deg F) 


V. PREDICTED PERFORMANCE CHARACTERISTICS 


In the actual application of the rotating heat pipe, various methods 
of cooling the condenser section are anticipated. The assumption is 
made here ,however, that the cooling mechanism employed is capable of 
setting the outside surface temperature of the condenser at a con- 
stant value, To Equation (38) is thus eliminated from the solution. 


Rearranging the remaining condensing limit equations gives 


3 
a 3 ee U4 1g 7 8334 
a SP (SM alg} eta omy" Gy 
a= 27 ky {hr lad + {Pot Fain $} eo 


These equations are a function of the physical properties of the working 
fluid, the thermal conductivity and the geometry of the pipe, and the 
operating conditions of rotational angular velocity and outside wall 
temperature TO0° 

A rotating heat pipe with the physical dimensions as shown in Fig. 
5 was chosen in order to demonstrate the predicted performance of a 
given pipe. The boiler section was modeled as shown in Fig. 5 in order 
to insure that an adequate supply of working fluid would be available 
for an actual pipe. Also with this model the rotating heat pipe boiler 


is closely approximated by the rotating boiler tested at Lewis Research 


32 


Center [7]. The results of the Lewis tests can be used to design 


and check the operation of the rotating heat pipe boiler. 


I ‘ Gh i G2 Gems 
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Figure 5) Simplified Model of a Rotating Heat Pipe 


The outside wall temperature of the condenser, TS was chosen as 

70 degrees Fahrenheit. For the saturation temperature of a parti- 
cular fluid, the fluid properties were determined. An IBM Dual 360 
digital computer was employed in iteratively solving the condensing 
limit equations (35) and (37). For a given angle, $, and rotational 
speed, W , the inside wall temperature, TO was varied until the 
heat rates of equations (35) and (37) agreed within 10 BTU per hour. 
It should be noted that the condenser length, 10 inches, and the 
radius of the boiler exit, 1.25 inches, were held constant while the 


angle @ was varied. 
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A, EFFECT OF FLUIDS 

The choice of the working fluid in a rotating heat pipe is dependent 
upon several physical properties of the fluid as well as the operating 
conditions desired. Various working fluids may be compared by the use 


of a dimensional parameter, N, defined as "the figure of merit." 


3 pe \/ 
y = {Antal * ith 
eal 
The heat rate will be proportional to the figure of merit of 
various fluids for a given saturation temperature. Table I gives the 


variation in N for three fluids at a saturation temperature of 122 


degrees Fahrenheit. 
TABLE I 


Figure of Merit 


Fluid N 
Water 20 
Ethyl Alcohol 5 
Freon 12 2.4 


It should be noted that for ordinary fluids, water is the most 
desirable working fluid. This is primarily due to the combination of 
a relatively high thermal conductivity (which is cubed) coupled with a 
high latent heat of vaporization. 

The vaporization curve of a given fluid will determine the operating 
temperature range of the fluid. The heat rate attainable will then be 
determined by the choice of fluid and the operating temperature of the 


pipe as shown in Fig. 6. It should be noted that in Fig. 6 the geometry 
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Figure 6 Heat Transfer Rate vs Saturation Temperature 
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of the rotating heat pipe is as shown in Fig. 5 with a constant rotation- 


al speed of 2700 RPM, and the half cone angle, 9, of 3 degrees. 


be EFFECT OF GEOMETRY 

For the given geometric configuration shown in Fig. 5 the half cone 
angle, $, was allowed to vary for a given operating condition and work- 
ing fluid. Figure 7 shows a plot of the heat rate versus the half cone 
angle for water at a saturation temperature of 212 degrees Fahrenheit. 
The rotational speed was allowed to vary from 300 to 3000 RPM. 

Figure 7 shows that a maximum heat rate is achieved at a particu- 
lar half cone angle for a given rotational speed (the dotted line). 
This peaking is due to the model used. The boiler exit diameter was held 
constant while the half cone angle, $, was varied. Since L.cos g, the 
condenser axial length was also held constant, the minimum radius of the 
condenser, R., varied. As $ increased R, decreased and thus the heat 
rate versus the half cone angle curve has a maximum point. 

Figure 7 illustrates that for the model chosen, as the rotational 
speed of the pipe is increased, the half cone angle necessary to main- 


tain the maximum heat rate decreases. 


C. EFFECT OF ROTATIONAL SPEED 

The rotational speed versus the heat rate is shown in Fig. 8 for 
various half cone angles, %. This curve is presented for water at 212 
degrees Fahrenheit in order to determine the most desirable operating 
speed range. One of the assumptions for the condensation on a truncated 
cone was that the centrifugal force is much greater than the gravity 
force (assumption h). For the model shown in Fig. 5, this condition is 


estimated to occur at approximately 700 RPM. 
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De OVERALL CAPABILITY 

A rotating heat pipe with the geometry shown in Fig. 5, operating 
at 2700 RPM, with a half cone angle of 1 degree, and using water as the 
working fluid was chosen to demonstrate the effects of the boiling, en- 
trainment, sonic, and condensing limits on the maximum heat rates at- 
tainable. 

The boiling limit curve was determined from equation (2), using 
both Kutateladze's value of K equal to .19 and Zuber's value of 7/24 
(.13), thus giving a range of heat rates for the boiling limit. Since 
the boiler surface of the rotating heat pipe may be approximated by a 
flat plate, for large radii of curvature, Zuber's predicted curve should 
be the most applicable for the model chosen. 

As mentioned earlier, no suitable theory is available to predict the 
entrainment limit, so that the experimental results of Hewitt, Lacey, 
and Nicholls have been used fae Their data was taken from a 1.25 inch 
diameter stationary tube in annular flow with a counter-current air- 
water contact length of 9 inches. They present flooding data for the 
vapor mass flow rate as a function of the liquid mass flow rate. The 
assumption of steady state operation in the rotating heat pipe allows 
the Hewitt, Lacey, and Nicholls data to be entered at the condition of 
equal liquid and vapor mass flow rates. The mass flow rate thus obtain- 
ed was then adjusted to agree with an average rotating heat pipe dia- 
meter of 2 inches and a counter-current steam-water flow condition. 

This adjusted mass flow rate of vapor was then used in equation (3) 


to give the entrainment limit. 


ig (3) 
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Equation (5) was used in computing the sonic limit. The assumption 
was made that choking occurred at the exit of the boiler section for a 
given saturation temperature in the condenser. 

The condensing limit was obtained by a solution of equations (35) 
and (37) as described on page 33. 

Figure 9 gives the plot of the rotational heat pipes heat rate as 
a function of the saturation temperature for the limits just discussed. 

For each of the limit curves shown in Fig. 9 there is an associated 
uncertainty. Even with this consideration, it is evident that the ro- 
tating heat pipe shown in Fig. 5 will be condenser limited, for a given 
saturation temperature. The upper limit on the heat rate attainable 
for a certain working fluid will Set) yep be governed by the ability 


of the container to withstand the saturation pressure. 


It is possible that a container could withstand a certain satura- 
tion pressure, and yet not be compatible with the associated wall 
temperature. 
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VI. CONCLUSIONS 


i. The film thickness and local heat transfer coefficients for 
condensation on a rotating, truncated cone may be determined by a 
Nusselt type of analysis. The resulting film thickness and local heat 
transfer coefficient are not constant as for the rotating cone and disc, 
but vary with axial position. 

2 The predicted performance for the heat pipe geometry chosen in 
this study, and operating with water as the working fluid, is con- 
denser limited. For a different geometry, the entrainment limit may 
become the dominant limit. The sonic and boiling limits should pre- 
sent no practical limitation on the rotating heat pipe. 

3% Of the ordinary fluids, water has the best characteristics 
for use as the working fluid. 

a The actual operation of the rotating heat pipe may be limited 
by the thermal conductance of the outside surface, (i.e., upon the 


method of cooling the rotating heat pipe). 
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VII. RECOMMENDATIONS FOR FURTHER STUDY 


A. ANALYTICAL 

ie An optimization study should be conducted to determine the 
effect of the internal geometry (boiler exit area and condenser shape 
as well as length) on the rotating heat pipe condensing limit. 

Ze The restrictions of negligible interfacial shear, negligible 
momentum changes in the condensate, and constant vapor pressure should 
be removed from the present analysis. This then would require a solu- 
tion of the continuity, momentum, and energy equations for the con- 
densate as well as the vapor. The equations of the condensate and the 
vapor would be coupled through the boundary condition that the shear at 
the interface, in each phase, must be equal. 

The swirl component of the flow of condensate should also be taken 
into account. This would then require that the problem be solved in 
three-dimensions, probably through the use of spherical co-ordinates. 

5% The exact location of the choked flow condition in the rota- 
ting heat pipe should be determined in order to more accurately chart 
the sonic limit. 

4. A study should be undertaken to determine accurately the en- 
trainment limit of the rotating heat pipe for the condenser geometry 


shown in Fig. 5. 


B. EXPERIMENTAL 
Lie An experimental model should be built and the condensing heat 
rate limit curve for the pipe should be experimentally verified for 


various working fluids. 


4,3 


2 A comparison of filmwise and dropwise condensation in the 
condenser section should be made. - 

St. A further extension of the experimental procedure should be 
to have a condenser wall with axial grooves in it. This then would 


test the contributions of interfacial shear and swirl on the condensing 


Limit. 
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APPENDIX A 


SOLUTION FOR THE FILM THICKNESS, & (X) 


Equation (30) may be written in the form 


z : dd 4. 
p, w h, sin d = 3{Ro+ Xsing} 5 au +25sindh a (A-1) 





Let ; 
Siu hy ginp 
me 
U=6 
Therefore 
8 G4 _K 
| ee | = A- 
RdU + {-5U--> fhe } aR=0 (A-2) 
Now let 
 @ £2 
a a he ain? 
This gives 
28085 90D. A-3 
aero ae (A-3) 
Integrating this gives 
Pose (A-4) 
- R 


fy 5 


where C is an arbitrary constant. Transforming equation (A-4) back 


into the original dimensions 


eg en (a5) 
3 SIN RotX sing 


pa 
eh) fo 


In order to determine the arbitrary constant, C, the boundary condi- 


tion 5 = 0 at X = 0 will be used. The constant, C, then is 


: a K 3/8 
C= Ro {- on “sing i (A-6) 


Now substituting the original value of K and the value of C into equa- 


tion (A-5), then rearranging, the final solution for 96 is 


o3 


ee ee wn 


pie w SIN $ he RotXsin 


The solution for condensation on a rotating cone, Sparrow and 
Hartnett's solution, can be easily obtained by the use of the boundary 
condition that 6 = 0 at X = 0 for Ro = 0. This condition would give, 
from equation (A-5), that C = 0. The solution for delta is then the 


solution obtained by Sparrow and Hartnett. 
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APPENDIX B 


SOLUTION FOR THE MEAN HEAT 


TRANSFER COEFFICIENT, hy, 


Equation (33) may be written in the form 


3 2 v4 
’ 2 {< k¢ Ae nia} (B-1) 
m L {2n+tsinp} bS os ' 3 
q4. 
it | Ror Xsin ie [" id. 
w/a x 
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+xXsinf 
Working with just the integral portion of equation (B-1) 
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Rot X sing 
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ein 
Equation (B-2) then becomes 
A cen, 5 
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Vi 
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It should be noted that 


Carrying out the integration therefore yields 


¢ I) 
ways “3 Ls ing} 





2 : —— r 8 —_ ve 
eer eae oe) 

a 

Q 

Finally 
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Substituting Z back into equation (B-l), the average heat transfer co- 


efficient, hw? is given by 
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